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ABSTRACT 

Aims. The principal aim of this project is to determine the structural parameters of the rapidly pulsating subdwarf B star PG 09 11 +456 
from asteroseismology. Our work forms part of an ongoing programme to constrain the internal characteristics of hot B subdwarfs 
with the long-term goal of differentiating between the various formation scenarios proposed for these objects. So far, a detailed 
asteroseismic interpretation has been carried out for 6 such pulsators, with apparent success. First comparisons with evolutionary 
theory look promising, however it is clear that more targets are needed for meaningful statistics to be derived. 

Methods. The observational pulsation periods of PG 0911-1-456 were extracted from rapid time-series photometry using standard 
Fourier analysis techniques. Supplemented by spectroscopic estimates of the star's mean atmospheric parameters, they were used as a 
basis for the "forward modelling" approach in asteroseismology. The latter culminates in the identification of one or more "optimal" 
models that can accurately reproduce the observed period spectrum. This naturally leads to an identification of the oscillations detected 
in terms of degree £ and radial order k, and infers the structural parameters of the target. 

Results. The high S/N low- and medium resolution spectroscopy obtained led to a refinement of the atmospheric parameters for PG 
0911+456, the derived values being 7^(1=31,940+220 K, logg=5.767±0.029, and logHe/H=-2.548±0.058. From the photometry it 
was possible to extract 7 independent pulsation periods in the 150-200 s range with amplitudes between 0.05 and 0.8 % of the star's 
mean brightness. There was no indication of fine frequency splitting over the ~68-day time baseline, suggesting a very slow rotation 
rate. An asteroseismic search of parameter space identified several models that matched the observed properties of PG 0911+456 
well, one of which was isolated as the "optimal" model on the basis of spectroscopic and mode identification considerations. All 
the observed pulsations are identified with low-order acoustic modes with degree indices f =0,1,2 and 4, and match the computed 
periods with a dispersion of only ~ 0.26 %, typical of the asteroseismological studies carried out to date for this type of star. The 
inferred structural parameters of PG 091 1+456 are 7^(1=31,940+220 K (from spectroscopy) , logg=5.777±0.002, M./Mo=0.39±0.01, 
logM,„,,/M.=-4.69±0.07, i?/i?o=0.133±0.001 and L/Lo= 16.4+0.8. We also derive the absolute magnitude Mv=4.82+0.04 and a 
distance <i=930.3+27.4 pc. 

Key words, asteroseismology, subdwarfs, pulsating stars 

1 . INTRODUCTION tion of their envelope mass near the tip of the first RGB, leav- 

^ , , ^ „ , , , . , , , ing them with insufficient fuel to ascend the asymptotic giant 

Subdwarf B (sdB) stars are evolved extreme honzontal branch ^^^^^^j^ ^^^g^ ^^^^^ ^^^^ j^^j^^^ exhaustion. Plausible forma- 

stars with atmospheric parameters in the 20,000 K < igfr < . , j n j • j » -i u lu TTI JtTvvtI 

Ar^r^^^r^ , ixT i i iTT^^^ ^. tion channcls were modelled in detail bv. e.g.. lHan et al.l (120021 

40000 K and 5.0 < log^ ^ 6.2 range dHebe^ta). They are ^ ,„d include binary evolution via a common envelope(cF) 

beheved to be composed of hehum-burnmg cores surrounded by j^^^^^ ^^^^j^ ^^^^^ j^^^ ^^^^.^^^ ^^^^^^ ^ 

thin hydrogen-rich envelopes and are characterised by ajiarrow ^^^-^^ ^^^^ ^^^^^^ ^^^^^ ^-^^^^^ evolutionary scenar- 

n^distnbution strongly peaked at ~ 48 ^DormanetalJ .^^ ^j^^^j^ j^^^^ ^ ^1^^^ ■ .^^ ^. 

119931) While It IS generally accepted tiiat they evolved from ^^^.^^ ^^^^^ evolution will produce sdB's in very 

the red giant branch (RGB) and con stitute the immediate pro- ^j^^^ ^. ^ ^^^^^^ j^^OF gives rise to much longer period bi- 

genitors of low-mass white dwarfs (j Bergeron et al. || l994 , the ^^.^^^ ^^.^^ ^^^^^ ^ ^^^^j^^ ^. ^ ^^^.^^^ 

details of their formation are not yet understood. It does how- ^j^^ ^^^.^ ^^^^ hydrogen-envelope thickness distri- 

ever seem clear that sdB progenitors lost a significant frac- ^^^.^^ Observational surveys focusing on radial velocity vari- 

* This smdy made extensive use of the computing facilities offered ^tions and the specti-oscopic detection of companions have re- 

by the Calcul en Midi-Pyrenees (CALMIP) project and the Centre centiy establis hed that at least half of the sd B population reside 

Informatique National de I'Enseignement Superieur (CINES), France, in binaries (e.g. lAUard et alJ(ll994 : IUlla & Thei ll ( 1998)), a sig- 

Some of the spectroscopic observations reported here were obtained at nificant fraction of them ha v ing short orbital peri ods from hours 

the MMT Observatory, a joint facility of the University of Arizona and to days jGreen et alj (ll997h : lMaxted et 311(120011) ). Accurate de- 
the Smithsonian Institution. 
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terminations of the stars' internal parameters on the other hand 
are harder to come by using traditional techniques; the mass has 
so far been measured only for th e very rare case of an eclips- 
ing binary dWood & Saffeiil 19991) . while the envelope thickness 
eludes direct study. 

Fortunately, a small fraction (~ 5 %) of sdB stars have 
been discovered to exhibit rapid, multi-periodic luminosity vari- 
ations on a time-scale of hundreds of seconds, thus open- 
ing up the possibility of using asteroseismology to constrain 
their internal parameters. S ince the near-sim ultaneous theoret- 
ical prediction (Charpinet et al. 19961 [19971) and observational 
discovery dKilkennv et al.l Il997l) of the so-called EC 14026 
stars, both the modelling and measurement of their pulsa- 
tional properties have come a long way (see [Fontaine et al.l 
( 1200 6^ for a review). Simulating the pulsation spectra of a 
large grid of sdB models in terms of low-degree, low-order 
p-modes and numerically determining the "optimal" model 
that best fits a series of observed periodicities has so far re- 
sulted in the asteroseismological determination of the inter- 
nal parameters for six EC 14026 pulsators: PG 0014+067 
drassard et al. 2 0^1), PG 10 47+003 (Charpinet et al. 2003), PG 
1219-1- 534 Charpinet et al.l "^OOSb), Feige 48 dCharpinet etH] 
I2005ah . EC 20117-4014 (.Randall et afl |2006a|), and PG 
1325+101 (ICharpinetet al.ll2006l) . These first asteroseismic re- 
sults show promising trends as far as matching the expected 
mass and hydrogen shell thickness distribution is concerned, 
however more targets are needed to start assessing the impor- 
tance and accuracy of the proposed formation channels. Here 
we present an asteroseismological analysis of the subdwarf B 
star PG 09 11 +456 based on photometry obtained with the new 
Mont4kccd at Mt. Bigelow, Arizona. In the next sections we de- 
scribe the observations and frequency analysis, followed by the 
asteroseismic exploitation of the target and a discussion of the 
internal parameters inferred. 

2. OBSERVATIONS 

2.1. Spectroscopy & atmospheric analysis 

PG 0911+456 was first identified as a hot star in the Palomar- 
Green survey (Green et al. 1986), and analysed in more detail by 
iKoen etal.l(ll999h using a medium-resolution (~0.8A/pix) spec- 
trum obtained on the 4.2-m William Herschel telescope on La 
Palma, Spain. Fitting a pure hydrogen line profile grid to the 
data suggested atmospheric parameters of Teff =3 1,900+200 K 
and log g=5. 82+0. 02 , while an independent analysis carried out 
by Rex Saff'er produced consistent values of refF=3 1,400 K and 
logg=5.80. This is typical for an EC 14026 star, and places PG 
0911+456 in the centr al part of the theoretical instability strip 
dCharpinet et al.ll200ll) . The authors also noted that the helium 
lines were weak for this object, and that there was no obvious 
spectroscopic signature of a companion. 

We obtained additional spectra with Steward Observatory's 
2.3-m telescope on Kitt Peak, Arizona, as well as the 6.5-m 
MMT as part of an ongoing program to determine the atmo- 
spheric characteristics of a large sample of both pulsating and 
non-pulsating sdB stars in a homogeneous way (for details see 
Green, Fontaine, & Chayer, in preparation). The Kitt Peak data 
cover the ~3600-5240 Arange at low wavelength resolution (~9 
A), while the MMT spectrum concentrates on the -4000-4950 
A region at a higher resolution of ~1 A. Both data sets were anal- 
ysed with the help of a detailed grid of non-LTE model atmo- 
spheres and synthetic spectra designed especially for subdwarf 
B stars. These were computed using the public codes TLUSTY 
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Fig. 1. Model fit (solid curve) to the hydrogen Balmer lines 
and helium lines available in our time averaged high S/N ratio, 
medium-resolution optical spectrum of PG 091 1+456. 



and SYNSPEC dHubenv & Lanzlll995l: iLanz & Hubenv"1995h. 
and in contrast to the models employed by Koen et al. (199^ 
incorporate helium but no metals. We show the best model 
fit to the available Balmer and neutral helium lines for the 
MMT and the 2.3-m spectra in Figures [T] and |2] respectively, 
and also indicate the atmospheric parameters inferred. Note that 
the values in brackets refer to the formal fitting errors. While 
both the eff'ective temperature and surface gravity derived from 
the low-resolution spectrum agree well with those reported by 
IKoen et a l] (fT999l) . the value of log g from the higher resolution 
data is significantly lower This is somewhat surprising, since in 
our experience the atmospheric parameters computed from both 
types of spectra are compatible as long as the S/N of the data is 
high enough. Examining the figures, we find the fits to both spec- 
tra to be satisfactory, however it is apparent that the MMT data 
are far noisier than those from Kitt Peak. Indeed, comparing the 
quality of the single medium-resolution spectrum obtained for 
PG 0911+456 to equivalent data f or previously analysed subd- 
warf B star s such as PG 1325+101 dCharpinet et al.ll2006l) or PG 
1219+534 dCharpinet et al.ll2005bl) reveals a significantly lower 
S/N (~78 compared to typically 200+) for the case of our tar- 
get. In contrast, our combined low-resolution spectrum of PG 
0911+456 is of outstanding quality (S/N~350). Given that the 
atmospheric parameters derived from the latter are also in ex- 
cellent agreement with those reported by Koen et al. ( 1999), we 
believe them to be more accurate than those obtained from the 
MMT spectra. This was also the conclusion reached by the ref- 
eree (U. Heber), who kindly re-analysed our Kitt Peak spec- 
trum for PG 091 1+456 using independent models that incorpo- 
rate metals at both solar and sub-solar abundances, and derived 
atmospheric parameters consistent with our low-resolution es- 
timates. He also noted that, in his analysis of MMT medium- 
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Fig. 2. Model fit (solid curve) to the hydrogen Balmer lines and 
helium lines available in our time averaged high S/N ratio, low- 
resolution optical spectrum of PG 091 1 +456. 



resolution spectra of about a dozen sdB stars, consistently lower 
surface gravities and temperatures were found than from their 
low-resolution counterparts, the latter being more reliable due 
to the wider wavelength coverage and higher number of Balmer 
lines sampled. 

Using the available estimates of the atmospheric param- 
eters for PG 0911+456 we compute a weighted mean of 
r^ ff=31,940±220 K and log ^=5.767+0.029. As already noted 
bv lKoen et al.l (Il999h . the atmospheric helium abundance is low 
even for a subdwarf B star at logHe/H=-2.548+0.058. 



a good sensitivity between ~3000 and 8000 A, with a quantum 
efficiency peak of about 95% at 4500 A. The read noise of the 
detector is 5.0 e" (at a readout rate of about 50 kHz per channel), 
the dark current is 16.6 e /pixel/hour at -130° C, the unbinned 
full well capacity is 131,000 e" and the device is operated at 
a gain of approximately 3.1 e"/DN. It is packaged in a Kovar 
tub and is steered with an Astronomical Research Cameras, Inc. 
Gen2 controller using two video channels. Target acquisition is 
performed with the help of the AzCam software developed by 
ITL and Steward Observatory. The Mont4kccd is equipped with 
a suite of filters (Bessell U, Harris BVR, Arizona I, as well as 
the broadband WFPC2 F555W and F814 and the Schott 8612 
filters, among othersEl). Mounted at the Cassegrain focus of the 
Mt. Bigelow telescope, the pixel scale is 0.14"/pixel, yielding a 
field of view of 9.7'x9.7'. 

The instrument was designed with differential photometry in 
mind, however the variety of filters available coupled with the 
excellent sensitivity of the chip and the large field of view make 
it ideal for many imaging projects. For the purpose of our obser- 
vations, the main challenge was to reduce the dead-time between 
exposures to an acceptable level. On-chip binning, normally 3x3 
pixels but occasionally 2x2, is ordinarily used to optimise the 
pixel size on the sky, improve readout speed, and save disk space. 
However, even 3x3 binning results in an overhead time of 18.2 
s, rather long when aiming to temporally resolve the multiple 
~150-200 s periodicities expected for an EC 14026 pulsator 
such as PG 0911+456. Our solution was to use 4x4 pixel bin- 
ning (0.6"/pixel) and read out only the middle 60% of the rows, 
reducing the overhead time to just under 10 s. For isolated stars 
we have in the past been able to obtain relative photometry accu- 
rate to 0.001-0.002 mags from aperture photometry of images 
with FWHM's as small as 1.75 pixels. Given that the seeing was 
never better than ~V during our observations of PG 091 1+456, 
this worked quite well for the observations presented here. The 
most difficult aspect of using a slightly degraded spatial resolu- 
tion with the Mont4kccd is maintaining the same effective fo- 
cus throughout the night, as the latter is mainly determined by 
instrumental flexure related to the pointing position rather than 
dome temperature. Since keeping the focus constant is neces- 
sary to avoid variations in the light curve, we developed a script 
to monitor the focus of the incoming images. 



2.2. Time-series pliotometry & frequency analysis 

PG 0911+456 was first identified as a member of the (then) 
newly disc o vered class of rapidly pulsating EC 14026 stars by 
lKoenetal.1 (Il999h . On the basis of -17 hours of Stromgen b 
band time-series photometry gathered at the 0.9-m McDonald 
Observatory on Mt. Locke, Texas, the authors were able to ex- 
tract 3 convincing periodicities between 150 s and 170 s (see 
Table 2 below) and estimate an approximate mean Stromgen b 
magnitude of ~14.6. 

Our follow-up photometry was obtained with the new 
Mont4kccd (also known as La Poune II) on Steward 
Observatory's 1.55-m Mt. Bigelow telescope in Arizona. Since 
this constitutes the first published science data obtained with 
the instrument, we give a brief overview of its characteris- 
tics and the observing strategy developed for fast time-series 
photometry. The 4kx4k CCD is a Fairchild CCD 486 de- 
tector with 4096x4097 15 micron pixels that has been pro- 
cessed for backside illumination at the Universtity of Arizona 
Imaging Technology Laboratory (ITL). Treatment with the ITL 
Chemiabsorption Charging and antireflection coating results in 



Table 1. Photometi-y obtained for PG 091 1+456 (2006-2007) 



Date (UT) 


Start time (UT) 


Length (h) 


20 Dec 


08:10 


05:07 


21 Dec 


07:40 


05:37 


15 Jan 


06:04 


07:21 


16 Jan 


05:30 


07:54 


17 Jan 


04:08 


09:19 


05 Feb 


02:55 


09:01 


10 Feb 


06:11 


04:05 


25 Feb 


02:37 


08:40 



We obtained a total of 57 hours of time-series photometry 
for PG 091 1+456 on 8 nights spread over a time baseline of 68 
days (see Table [1] for details). Owing to the relative faintness 
of the target, the exposure time was set to 30 s throughout the 
campaign, and we employed the wide-band Schott 8612 filter 



' The Mont4k;ccd QE curve as well as the filter trans- 
mission curves can be found on the following webpage: 
http : //j ames . as . arizona.edu~psmith/6 1 inch/instruments .html , 
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Fig. 3. All light curves obtained for PG 0911+456 on the Mt. 
Bigelow 1 .55-m telescope. The data have been shifted arbitrarily 
along the x and y axes for visualisation purposes. From top to 
bottom, the curves refer to the nights of 20 Dec, 21 Dec (2006), 
15 Jan, 16 Jan, 17 Jan, 05 Feb, 10 Feb and 25 Feb (2007). For 
details see Table 1 . 



to maximise throughput. The data were reduced using standard 
IRAF reduction aperture photometry routines, except that we set 
the photometric aperture size separately for each frame to 2.5 
times the FWHM in that image. We computed differential light 
curves of PG 0911+456 on the basis of 6 suitable comparison 
stars distributed as symmetrically as possible around the target. 
Figure [3] shows all the nightly light curves obtained; the beat- 
ing between different modes is clearly visible. However, the re- 
markable quality of the photometry is better appreciated from 
the expanded view of a single night's light curve as displayed in 
FigurelH The relative noise level is reminiscent of that achieved 
on much larger telescopes, such as the 3.6-m CFHT or the 3.5- 
m NTT (see, e.g. Figure 2 of Brassard et al., (2001), Figure 4 



of ICharpinet et al.1 (l2005bl) and Figure 2 of iBilleres & Fontaing 
(Hooi) when taking into account the differences in target bright- 
ness. We attribute this partly to the excellent sensitivity of the 
CCD and our optimised data pipeline, but even more importantly 
to the relatively high S/N achieved in each exposure and to the 
excellent spatial distribution of comparison stars around the tar- 
get. Note that a high data quality can be achieved under less 
than perfect atmospheric conditions including transparency vari- 
ations, thin cirrus and a bright Moon, rendering this technique 
very robust. 

The time-series photometry gathered for PG 091 1+456 was 
analysed in a standard way using a combination of Fourier anal- 
ysis, least-squa re fits to the light cur ve, and pre-whitening tech- 
niques (see e.g. IBilleres et al. ( l2000l) for more details). Figure|5] 
shows the Fourier spectra of the individual nightly light curves. 



B 



00 ?'.^/.^X^^^'WV\^V-*,-«*V.^^^^^rt^V*dA^V 

-0.05 L . : : : : 1 : : : ' : L 



0.05 — ' — I — — 

0.00 

0.05 L — : — I — : — : : — I — ^ 

14000 15000 



0.05 
0.00 
-0.05 



0.05 
0.00 

-0.05 



37000 SBOOO S9000 30000 31000 33000 33000 



Time (seconds) 

Fig. 4. Expanded view of the longest light curve obtained on 16 
Jan 2007 in units of fractional brightness intensity and seconds. 
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Fig. 5. Montage of the nightly Fourier transforms, shifted arbi- 
trarily along the y-axis for visualisation purposes. 



which all clearly show a dominant peak at ~ 6.4 mHz as well as 
a few lower amplitude frequencies in the 6-7 mHz range and an 
apparently isolated weak oscillation at ~ 5.2 mHz. The power 
spectrum beyond 13 mHz out to the Nyquist frequency of 26 
mHz is entirely consistent with noise and hence not illustrated. 
Given the poor duty cycle of the data over the long time baseline 
(3.5 %), it was not obvious that combining all the light curves to 
give a total Fourier transform would be fruitful and produce reli- 
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Fig. 6. Fourier transform of the entire data set in the 5-7 mHz 
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by the dominant, the two strongest and finally all 7 frequencies 
detected above 4cr. 



able peaks (rather than windowing artefacts). We therefore per- 
formed a Fourier analysis for both the full dataset and a subset 
of three consecutive nights (14-16 January) and compared the 
results. Somewhat gratifyingly, the frequencies and amplitudes 
extracted were the same within the measurement uncertainties 
and the noise level decreased slightly when adding extra data 
sets. This has not always been the case for past observational 
campaigns of EC 14026 stars, where combining light curves 
well separated in time or of inhomogeneous quality signifi- 
cantly degraded the pulsati on spectrum or pro d uced incoherent 
osci llation peaks (see, e.g. Stob ie et al] (Il997h : ICha rpinet et al.' 
( l200 5b a)). The successful merging of different light curves for 
PG 0911+456 was most likely possible due to the high S/N of 
the individual light curves as well as a relative homogeneity be- 
tween the data sets. In what follows, we will focus on the results 
obtained from all the data, since they correspond to the highest 
frequency resolution. 

The Fourier transform of the combined light curve is dis- 
played at the top of Figure |6] followed by successively pre- 
whitened spectra in the lower part of the plot. In each case, the 
number of periodicities extracted is indicated above the spec- 
trum. We detected a total of 7 oscillations down to a threshold of 
4 times the noise level in the frequency range plotted (0.07 %). 
The quality of the pre-whitening can be assessed from Figure 
|7] which shows the original Fourier transform together with a 
noise-less spectrum reconstructed on the basis of the 7 extracted 
periodicities and the point-by-point difference between the two. 
While the latter indicates the possible presence of additional, 
unidentified frequencies, we feel that we have reached the de- 
tection limit for convincing oscillations. The periods, frequen- 
cies, and amplitudes of the pulsations extracted from the data 
are listed in Table |2] It is striking how well the three domi- 
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Fig. 7. Comparison between the Fourier transform of the entire 
data set (Full spectrum), and that reconstructed on the basis of 
the 7 periods, amplitudes and phases derived from least-squares 
fitting to the light curves (Model spectrum). The lower curve 
shows the point-by-point difference between the actual FT and 
the model FT. 

Table 2. Oscillations extracted from the combined light curve of 
PG 0911+456. The uncertainty on the frequencies is estimated 
to be about 1/lOth of the frequency resolution, i.e. around 0.02 
fiHz, while the error on the amplitude is derived from the least- 
squares fit to the light curve. The frequencies and amplitudes 
detected bv lKoen et al.l d 19991) are also listed. 



Rank Period 


Frequency 


Amplitude 


./koen -^Koen 


(s) 


0"Hz) 


(%) 


(y"Hz) (%) 



4 


192.551 


5193.42 


0.082±0.005 






6 


168.784 


5924.73 


0.054+0.005 






2 


165.687 


6035.49 


0.223+0.005 


6036.2 


0.44 


3 


161.554 


6189.89 


0.098+0.005 


6190.9 


0.24 


5 


157.581 


6345.92 


0.080+0.005 






1 


155.767 


6419.85 


0.755+0.005 


6419.3 


0.61 


7 


149.027 


6710.20 


0.050+0.005 







nant peaks correspond to the oscillations detected bv'Koe n et al.l 
( 1999) (/Koen) some 9 years previously. Even the amplitude rank 
order has remained the same, although the two weaker pulsations 
seem to have diminished somewhat with respect to the main peak 
(note that the absolute values of the amplitudes cannot be read- 
ily compare d since the bandpass and frequency resolution of the 
iKoen et all ( [1999) data are different to our own). Such ampli- 
tude variations are by no means uncommon in EC 14026 stars, 
and can be readily interpreted either in terms of the beating of 
unresolved modes (e.g. very closely spaced components of a ro- 
tationally split multiplet) or significant changes in the intrinsic 
amplitudes of the modes. In this case, we lean towards the lat- 
ter explanation since it is likely that at least one of the dominant 
peaks is associated with a radial mode and therefore not subject 
to unresolved rotational splitting. And indeed, the two highest 
amplitude pulsations are identified with f-Q modes in the aster- 
oseismological analysis presented below. 

The absence of frequency splitting over the campaign base- 
line of over two months implies that PG 0911+456 is a slow 
rotator with a rotation period Prot ^ 68 d. It is of course also 
possible that the components of any multiplets have amplitudes 
below the detection limit, however this seems unlikely given the 
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low noise level of the data. The inferred long rotation period cou- 
pled with no sign of a companion being detected from either the 
time-averaged spectroscopy or the 2MASS J-H colour suggests 
that our target is a single star This hypothesis is supported by 
radial velocity measurements obtained at Steward Observatory's 
2.3-m Bok telescope on Kitt Peak as part of an ongoing spec- 
troscopic survey of sdB stars (Green et al, in preparation). No 
detectable velocity variation was found from ten spectra of PG 
09 1 1 +456 taken at irregular intervals over a period of three years 
(2002-2005). While at least half of all subdwarf B stai's appear 
to reside in binary systems, these are thought to mostly have 
short orbital p eriods (hours to days) and white dwarf secondaries 
dGreen etal.l ll997; M axte d et al. 2001) or rather long periods 

1 hundreds to thous ands of days) and main sequence companions 
Saffer et ani200lb . Assuming tidally locked rotation, the low 
noise level of our photometry should reveal rotational frequency 
splitting in the former case, and a spectroscopic signature of a 
cooler companion in the latter Depending on the inclination of 
the system, the binary motion of the subdwarf may also manifest 
itself in a Doppler shift of the pulsation frequencies on the time- 
scale of the orbital p eriod (as detected for t he long-period binary 
EC 201 17-4014, see lRandall et all (l2006ah), of which a co mpar- 
ison of our data with those obtained by Koen et aT] (Il999i) show 
no indication. 



3. Asteroseismic Analysis 

Our asteroseismological study of PG 09 1 1 +456 was carried out 
using the well-known f orward method, describe d in detail by e.g. 
iBrassard et"an ( 1200 lb : ICharpinet et all (l2005b). It consists of a 
double-optimisation procedure that first determines and quanti- 
fies the best match between the set of observed periodicities and 
those calculated for a given model, and subsequently searches 
parameter space for the model (or models) that can reproduce 
the data most accurately. The solution is referred to as the "op- 
timal model", and corresponds to the region of parameter space 
where the goodness-of-fit merit function S ^ reaches a minimum. 
Note that 5 ^ is not a standard since no weighting is applied, 
and is given by 



(1) 



where P'^^^^ is one of the n periodicities observed and ^'J;,^^ is the 
theoretical period that matches it best. 

Over the years we have developed various sets of numerical 
tools that allow us to apply this scheme to isolate best period 
fitting models of subdwarf B stars for asteroseismology. Two 
independent packages, fortunately leading to the same results, 
are now operational, one in Montreal and the other in Toulouse. 
Both are built on the simple requirement of global optimisa- 
tion, i.e. models aiming to provide a good asteroseismological 
fit must be able to match all the observed periods simultane- 
ously. Mode identification is normally not assumed a priori, but 
instead appears naturally as the solution of the global fitting pro- 
cedure. This provides a convenient cross-check of the "optimal" 
model's robustness in the case where independent constraints on 
mode identification are available, e.g. from rotational splittin g 
or multi-colour photometry (see e.g. ICharpinet et al.l (l2005ah . 
Brassard et al. (2007) in preparation). 

The theoretical period spectra are computed on the basis of 
our so-called "second-generation" models (see Charpin et et al.l 
(119971 12001h ). These models are static envelope structures in- 
corporating microscopic diffusion under the assumption of an 



equilibrium having been established between gravitational set- 
tling and radiative levitation. Since it is diffusion that creates the 
iron reservoir responsible - through the /c-mechanism - for pul- 
sational instabilities in sdB stars, its detailed treatment is essen- 
tial when predicting the excitation of modes. Moreover, micro- 
scopic diffusion changes the stellar structure sufficiently to mod- 
ify the pulsation periods themselves, and is therefore a neces- 
sary ingredient even for the calculation of adiabatic periods. The 
internal structure of the models is specified by four fundamen- 
tal input parameters: the effective temperature Teff, the surface 
gravity log g, the total stellar mass M, and the fractional thick- 
ness of the hydrogen-rich envelope logq'(H) = log [M(H)/M,]. 
The latter parameter is intimately related to the more commonly 
used quantity Mf,,,,,, which corresponds to the total mass of the 
H-rich envelop^ Both adiabatic and non-adiabatic oscillation 
modes are computed for each model using two efficient codes 
based on finite element techniques. The first solves the four adi- 
abatic oscillation equations and constitutes an intermediate but 
necessary step to derive estimates for the periods and eigenfunc- 
tions that are then used as first guesses for the solution of the 
full, non-adiabatic set of oscillation equations. The latter com- 
putes a number of quantities for each mode with given degree 
€ and radial order k, most importantly those that can be directly 
compared with observed quantities, such as the period f ,/, and 
the stability coefficient cr/. If cr/ is negative, the mode is stable, 
while if it is positive the mode is excited, and may therefore be 
observable if its amplitude is large enough. This can be used as a 
cross-check for the reliability of the optimal model inferred: ob- 
viously, all the periodicities detected should be associated with 
unstable theoretical modes. Other quantities output by the non- 
adiabatic code include the kinetic energy E and the rotation co- 
efficient Cki (see Table |4]i. The former is of interest mainly for 
theoretical studies, e.g. when investigating mode trapping, and 
the latter becomes important only when rotational splitting is de- 
tected; therefore they are of secondary interest for our purposes. 

The observations of PG 0911+456 provide us with 7 har- 
monic oscillations to compare with theoretical predictions. This 
seems a reasonable number for attempting an asteroseismologi- 
cal analysis, since similar studies have been carried out on the 
basis o f anywhere between 3 ( EC 201 17-4014. [Randall et al.l 
([2006ah ) and 13 (PG 0014+067. iBrassard et all (120011) ^ period- 
icities. Nevertheless, the model finding algorithm uncovers sev- 
eral families of models that can account for the detected pe- 
riod spectrum similarly well, complicating the identification of a 
convincing "optimal" model. By carefully exploring parameter 
space under different assumptions as to the allowed degree in- 
dices of the modes observed we nevertheless isolate a promising 
solution. 



3.1. Search for the optimal model for PG 091 1+456 

We started the search for the optimal model by keeping any a 
priory assumptions to a minimum. As usual, the pulsation cal- 
culations were effectuated for all modes with degree indices 
< ^ < 4 in a period range amply covering the observed oscilla- 
tions (100-300 s). Higher degree modes are thought to have am- 



- Note that the parameter M„„, commonly used in extreme horizon- 
tal branch stellar evolution includes the mass of hydrogen contained in 
the thin He/H transition zone, whereas the parameter M{H) used in our 
envelope models does not. They can be related with log [Mj,„,/A^*] = 
log [M{H)IMt] + C, where C is a small positive term slightly depen- 
dent on the model parameters that can be computed from the converged 
model using the mass of hydrogen present in the transition zone itself. 
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Fig. 8. Left Panel: slice of the S ^ function (in logarithmic units) along the log g - T^s plane at fixed parameters M» and log q(li) set to 
their optimal values found for the best-fit solution Model 1 (M»=0.40 Mq and logq'(H)=-4.11). The solid line rectangle represents 
our mean spectroscopic estimate of the atmospheric parameters with their uncertainties. Also indicated are the 1, 2 and 3 cr limits 
on the minimum (white dashed contours). The errors on the asteroseismologically derived parameters are estimated from the 
semi-axes of the 1 cr contour. Right Panel: slice of the S ^ function (in logarithmic units) along the M, - log ^(H) plane at fixed 
parameters \ogg and Teff set to their optimal values found for the best-fit solution Model 1 (logg=5.686 Mq and reft=31,940 K). 
The 1, 2 and 3 cr limits on the minimum are once again indicated. 



plitudes too small for them to be detectable due to cancellation 
effects when integrating over the visible disk of the star Indeed, 
the observability of modes with degrees higher than 2 was long 
debated, but the period density and distribution uncovered in cer- 
tain EC 14026 star s (e.g. KPD 193 0+2752, lBill6res et al.i (l2000h 
or PG 1325-HlOl. fCharpinet etall (l2006h ) forced the inclusion 
of higher degree modes. In the case of PG 09 1 1 -1-456 this is not 
strictly necessary as the pulsation periods observed can be ac- 
counted for solely in terms of modes with / < 2 (see below). 
However, we feel that it would be counterproductive to exclude 
potentially interesting models from the outset. 

Throughout the exploration of parameter space, we fixed the 
effective temperature to the mean spectroscopic value of 31,940 
K, thereby decreasing the number of free parameters from 4 to 3. 
In our experience, this facilitates and speeds up the identification 
of potential "optimal" models. Given the weak temperature de- 
pendence of the pulsation periods, Tetj is in any case determined 
more accurately on the basis of spectroscopy than asteroseismol- 
ogy (the inverse holds true for log g, which can be constrained 
far more tightly from asteroseismology by virtue of its strong im- 
pact on the period spectrum). And indeed, test calculations car- 
ried out at fixed temperatures at the limits of the spectroscopic 
uncertainties (employing Teff =32,200 K and Teff =3 1 ,800 K) pro- 
duced the same families of optimal models, albeit with a slight 
(~ 0.01 Mq) shift in mass in each case. The compensation of an 
increase in temperature with an increase in mass is well under- 
stood in terms of parameter degeneracy and has been observed 
in the asteroseismi c exploration of most ot her EC 14026 stars 
(e.g. PG 1 2 1 9 -h5 34. ICharpinet et all (12005b ')) . It implies that the 
two parameters cannot be optimised independently, obliging us 
to fix or at least constrain the temperature if the mass is to be 
determined to any accuracy. As for the remaining model input 



parameters, we sandwich the value suggested from spectroscopy 
for the surface gravity, searching the range 5.50 < logg < 5.80, 
while possible values of the stellar mass and the hydrogen en- 
velope thickness are constrained from evolutionary models to 
0.30 < M,/Mo < 0.7 and -5.20 < logq{H) < -2.00 (see 
ICharpinet et al.l (l2005bb for more details). 

Within the search domain specified, our algorithm identified 
several families of models that present good period matches. The 
majority of these were discarded upon closer inspection, either 
because of obvious inconsistencies with the spectroscopic values 
of log g, or because the mode identification inferred was deemed 
inappropriate (e.g. the dominant pulsation was assigned {=3). 
However, two of the best-fit models (i.e. those that matched the 
observed periods with the lowest values of S ^) retained our at- 
tention as potential "optimal" models and are discussed in more 
detail in what follows. 

The first minimum in the goodness-of-fit merit function 
(at 5^ ~ 0.63) corresponds to a model with Teff - 31,940 
K (fixed from spectroscopy), log ^=5.686, M»=0.40 Mq and 
logq(H)--4.ll (hereafter referred to as Model 1). Figure|8]il- 
lustrates the complex behaviour of the S ^ hyper-surface in the 
vicinity of this solution: the left hand panel shows a slice along 
the log g - Teff plane with the mass and envelope thickness set 
to their optimum values, whereas the right hand panel displays a 
cut through M, - log ^(H) space with the temperature and surface 
gravity fixed to the Model 1 solution. As indicated on the scale 
to the right of the main plots, areas shaded in dark blue indicate a 
relatively good match between the observed and theoretical pe- 
riods, while those coloured in red represent comparatively bad 
fits. From the left hand panel we can see that the good-fit mod- 
els are concentrated in two elongated S ^ "valleys", one centred 
around logg ~ 5.68 (Model 1), and the other around logg ~ 
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Table 3. Mode identification inferred from Model 1 and Model 
2. 
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5.61. The second of these is quite clearly inconsistent with the 
mean spectroscopic estimate (indicated by the black rectangle) 
and can be rejected immediately. Model 1 on the other hand may 
lie well outside the spectroscopic box shown, however a closer 
look reveals that the surface gravity inferred is quite close to that 
derived from our MMT spectrum. Despite our conclusion that 
the latter is less accurate than the higher log g value suggested 
by both the Kitt Peak data and K oen et a l. ( 1999) (see Section 
2.1), it seems foolhardy to reject this solution immediately. 

The second possible family of optimal models (with ~ 
1.64) is centred at reff=31,940 K (fixed from spectroscopy), 
logg=5.777, M.=0.39 and log^(H)=-4.78 (hereafter referred 
to as Model 2). Although the value of is more than a factor 
of 2 higher than that for Model 1, the period match is neverthe- 
less perfectly acceptable at < AP/P >~ 0.26 % or, equivalently, 
an absolute dispersion of 0.41 s, and quite typical of the preci- 
sion achieved in asteroseismological studies to date. Moreover, 
the value of the surface gravity determined is nearly identical to 
our most reliable spectroscopic estimate. Seeing as though both 
Model 1 and Model 2 predict all the theoretical modes matched 
with observed pulsations to be unstable, non-adiabatic theory 
does not further discrimination between the two. Instead, we turn 
to the mode identification inferred. 

Table [3] lists the observed pulsations in order of decreasing 
amplitude rank, together with the resulting mode identification 
in terms of degree { and radial order k for the two solutions. 
In the case of Model 1, the degree indices assigned to the ob- 
served oscillations are not entirely convincing. While the domi- 
nant pulsation is associated with a radial mode as would be ex- 
pected, both on the basis of visibility arguments and partial mode 
identification from multi-colour photometry for other EC 14026 
stars (the main pulsations in both KPD 210 9-H4401 and Balloon 
090100 001 were s hown to be radial, see Randall et all (l2005b 
and .Fontaine et al.l (2006J respectively), the amplitude rank 3 os- 
cillation is identified with an ^=3 mode. Detailed model atmo- 
sphere calculations have recently shown that such modes have 
negligible amplitudes an d are in fact geometrically less visi- 
ble than those with E-A jRandall et al.l (120051) : see also Figure 
|9|l. Moreover, multi-colour photometry of pulsating subdwarf B 
stars has never indicated the presence of ^=3 modes at detectable 
amplitudes, and specifically excluded this possibility f or some 
27 oscillations observed in 7 different stars (see Trembla v et al.l 
(120061) . iJefferv et al.l (l2005h ). Note that these pulsations show a 
very different amplitude-wavelength behaviour to modes with 
^=0, 1,2,4 and should therefore be readily identifiable if they 
are found. Modes with {-A on the other hand have been in- 
ferred in a few cases, e.g. for KPD 2109-H4401, HS 0039H-4302 
( iJefferv et al.1 i2004 and Balloon 090100001 (Brassard et al. 
2007, in preparation). 



j=3 _ [_ ' 
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Fig. 9. Variation of the disk-integrated pulsational amplitudes 
with inclination angle / for m-Q modes with i-0 (continuous 
line), {=1 (dot-dashed), £=2 (dotted), ^=3 (dashed - hardly visi- 
ble due to the extremely low amplitudes) and f-4 (long-dashed) 
for an atmospheric model of PG 0911+456. The amplitudes 
were normalised with respect to the radial mode. Note that the 
intrinsic amplitudes were assumed to be the same for all modes. 



Model 2 seems rather more convincing in terms of the de- 
gree and radial orders inferred. For one, the amplitude hierarchy 
goes strictly with £: the dominant two modes are found to be 
radial, followed by two quadrupole modes, two pulsations with 
f-4 and, finally, a lone dipole mode. Secondly, with the excep- 
tion of the weakest oscillation, the modes observed for each de- 
gree correspond to the fundamental mode and the first overtone, 
indicating bands of instability rather than isolated frequencies 
excited to observable amplitudes. Although our pulsation calcu- 
lations rely on linear theory and consequently cannot predict the 
amplitudes of modes, it is more intuitive to imagine that within a 
certain frequency range all modes should be excited to a compa- 
rable intrinsic amplitude. And, indeed, asteroseismological stud- 
ies of other EC 14026 pulsators find that the majority of observed 
modes, particularly those with relatively high amplitudes, lie at 
the low-frequency end of the p-mode branch. 

One potential shortcoming of both Model 1 and Model 2 is 
the conspicuous absence of dipole modes among the higher am- 
plitude oscillations listed in Table[3] According to canonical wis- 
dom, the amplitude at which a mode is observed should decrease 
rapidly with increasing degree index due to cancellation effects 
when integrating over the visible disk of the star. However, de- 
tailed calculations employing full model atmospheres (taking 
into account limb darkening, as well as radial, surface gravity 
and temperature variations over a pulsation cycle) show that this 
i s not necessarily true. It was already shown by Randall et al] 
(|2005|) that {-4 modes are more visible than those with £-3. 
Here, we go one step further and investigate the relative am- 
plitude of modes with different degree indices as a function of 
the inclination angle / (defined as the angle between the line- 
of-sight and the pulsation axis) employing the modified model 
atmosphere code for subdwarf B stars developed in Montreal. 
This code incorporates several programs utilising the same bank 
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Fig, 10. Left Panel: slice of the S ^ function (in logarithmic units) along the log g - Teff plane at fixed parameters and log q(H) 
set to their optimal values found for the best-fit solution (M^-0.39 M© and log q(H)--4.1S). The solid line rectangle represents our 
mean spectroscopic estimate of the atmospheric parameters with their uncertainties. Also indicated are the 1, 2 and 3 cr limits on the 
minimum (white dashed contours). The errors on the asteroseismologically derived parameters are estimated from the semi-axes 
of the 1 cr contour. Right Panel: slice of the function zoomed in on the minimum. 
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Fig. 11. Left Panel: slice of the function (in logarithmic units) along the M, - log qiH) plane at fixed parameters logg and Teff 
set to their optimal values found for the best-fit solution (log^=5.777 and retf=31,940 K). Also indicated are the 1, 2 and 3 cr 
limits on the S' minimum (white dashed contours). Right Panel: slice of the function zoomed in on the minimum. 



of full m odel atmospheres an d was used and described exten- 
sively by iRandall et al.l (l2005h . Given values for logg and Teff 
as well as an inchnation angle the wavelength bandpass of 
interest, and the characteristics of the oscillation mode, it out- 
puts the emergent perturbed Eddington flux at user-specified 
points in time, and consequently enables the determination of 



relative pulsational amplitudes. In our calculations, we fixed 
the atmospheric parameters to those derived for PG 0911-1-456 
from spectroscopy, set the bandpass to white light received at 
a representative observing site (in this case Kitt Peak National 
Observatory), and calculated the integrated amplitudes expected 
for low-degree {0 < { < 4) modes for inclination angles from 
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to 90° assuming they have the same intrinsic amplitudes. Note 
that the azimuthal order m was set to as appropriate for a slow 
rotator. The results are illustrated in Figure |9] which shows the 
observable amplitudes expected for i-l,2,3 and 4 modes nor- 
malised relative to the radial mode (the amplitude of which is 
of course independent of the inclination) as a function of incli- 
nation /. It is quite striking that the canonical f=0,l,2,... ampli- 
tude hierarchy does not always hold true: in particular, dipole 
modes experience a rapid decrease in visibility with increasing 
inclination, and quadrupole modes have a higher integrated am- 
plitude above -80°, while being virtually undetectable at ~55°. 
As expected, £=3 modes have very low visibility regardless of 
the inclination angle. 

Returning to the mode identification inferred from our two 
models, we find that the absence of high-amplitude dipole modes 
could be explained by a very large inclination angle (/ > 80°). 
Moreover, the amplitude-degree index hierarchy inferred for 
Model 2 then fits in very nicely with the plot, that for Model 
1 less so due to the high-amplitude £-3 mode identified. Of 
course, we are working under the assumption that the intrinsic 
oscillation amplitudes of the modes are similar, which in the ab- 
sence of non-linear pulsation theory we have no means of ver- 
ifying. Nevertheless, we believe Figure |9] carries the important 
message that the canonical simplified {-0,1,2... amplitude hier- 
archy is not always valid, and that the visibility of a mode is 
strongly dependent on the inclination at which it is observed. 

Having established that modes with f=4 are more likely to be 
visible than those with {=3, we conducted a second search for an 
optimal model specifically excluding the latter. The ranges con- 
sidered for the input parameters were kept the same as during 
the first search. Perhaps unsurprisingly considering the quality 
of the fit, the global minimum was found to be identical with 
the Model 2 isolated during the first sweep of parameter space. 
Fi gures [TOl andfTTI show the corresponding slices of the func- 
tion along the log g - Teff and the - log ^(H) plane respec- 
tively, with the remaining parameters fixed to their optimal val- 
ues. As already pointed out above, the value of log^ inferred 
is in excellent agreement with the mean spectroscopic estimate. 
It is also apparent that there are two shallow valleys along 
the effective temperature axis, one centered on Model 2, and the 
other around \ogg ~ 5.66, close to Model 1. This is not a co- 
incidence; in fact, we demonstrate below that Model 1 and 2 
represent analogous solutions, with a mode identification shifted 
by k--l from the lower to the higher gravity valley. Note 
that in an extension of this "mode jumping" phenomenon the 
third minimum at \ogg ~ 5.61 in Figure [8] is associated with 
a ^=-1-1 shift from Model 1. This behaviour can be understood 
in terms of a relative uniformity of the period distribution for a 
given { value and the strong period dependence on log g. When 
the surface gravity of a model is slowly moved away from that 
of the optimal model, the theoretical period spectrum shifts com- 
pared to that observed, and the quality of the fit degrades at first. 
Once the change in log g becomes large enough, the S ^ function 
slowly decreases, converging towards the next solution, with a 
mode identification shifted by 1 in radial order Luckily for as- 
teroseismology, the low-order /7-modes detected in EC 14026 
stars are quite sensitive to the internal composition of the star, 
and the modes are normally not spaced too uniformly in k, en- 
abling a discrimination between different solutions on the basis 
of 5 ^ in conjunction with the spectroscopic constraints on log g 
and Tsff (see Figure[T2li. In the high-order asymptotic regime the 
story is very different, and this is in fact one of the main prob- 
lems when quantitatively i nterpre ting slowly pulsating subdwarf 
B stars (see lRandall et alj (l2006bj) ). 



In order to assess the significance and relative importance 
of Model 2, we carried out several additional searches in 3- 
dimensional parameter space under artificially imposed con- 
straints. First, we limited the acceptable degree indices to ^ < 
2 to determine whether the observed period distribution could 
be recreated solely in terms of the traditionally expected modes. 
While we found a physically viable model, the quality of the fit 
(S^ ~ 38) was so severely degraded compared to the previous 
models that it could not be seriously considered as an optimal 
model. Next, we forced the degree index identification inferred 
from Model 1 and 2 onto the higher and lower gravity valleys 
respectively, and were able to recover local minima very close to 
the other model, albeit at a much degraded quality-of-fit. In both 
cases, the mode identification was shifted in k. Model 2 start- 
ing at k-Q and Model 1 at A:=l in a convincing illustration of 
"mode jumping". Finally, we imposed the canonically expected 
amplitude -degree hierarchy on mode identification during the 
optimal model search, limiting the higher amplitude modes to 
f <2, while additionally allowing {-4 for the two lowest ampli- 
tude modes. From this we found a model very similar to Model 
2, at a degraded but still acceptable quality-of-fit (S^ ~ 9.8). We 
believe that this confirms the fundamental robustness of Model 
2, and thus adopt it as the optimal model for PG 091 1+456. 

3.2. Period fit for the optimal model 

The optimal model isolated for PG 09 1 1 +456 provides an excel- 
lent match to the 7 independent periods clearly identified in this 
star Details on the period fit and mode identification are given in 
Tableland are graphically represented in Figure [121 As briefly 
mentioned earlier, the quality of the period fit at AP/P ~ 0.26 
% is typical of that achieved for asteroseismological analyses of 
other EC 14026 stars; the dispersion is however still an order of 
magnitude higher than the accuracy to which the observed peri- 
ods can be measured. This has been the case for most EC 14026 
stars analysed to date and is attributed to remaining imperfec- 
tions in the equilibrium models. It is of course one of the goals of 
our asteroseismic studies to fine-tune these models, a quest that 
we are currently actively working on. Nevertheless, we believe 
that the fact that we were able to isolate an optimal model that a) 
is in accordance with spectroscopy, b) can accurately reproduce 
the observed period spectrum and c) predicts all observed modes 
to be unstable confirms the basic validity of our sdB star models. 

From Figure [12] it is apparent that the observed periods 
are exclusively identified with modes clustered at the high- 
period end of the theoretical p-mode branch. This behaviour 
is reminiscent of that uncovered for other EC 14026 pulsators 
such as PG 0014+067 (^ssard et al.l l2d01) and PG 1219+534 
dCharpinet et al.l l2005bh . and seems to support the idea that 
energy powering up the pulsations may be preferentially dis- 
tributed among the modes of low radial order. These modes 
would then reach observable amplitudes more easily (depending 
of course on the geometric mode visibility) and would therefore 
be detected in the light curves more readily. It should be noted 
however, that some E C 14026 pulsators (e.g. PG 1325+101, 
ICharpinet et al.l ([2006')) also show isolated modes at lower pe- 
riods, with apparent gaps in the radial order distribution. These 
cannot be fully understood without non-linear pulsation theory, 
but it seems probable that some of the "missing" modes will 
eventually be detected. In the past, observations at higher S/N 
have always revealed extra oscillations compared to noisier data, 
supporting the hypothesis that sdB stars often excite many more 
modes than are uncovered, albeit at lower observable ampli- 
tudes. Apart from geometric visibility, the observability of oscil- 
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Table 4. Period fit for the optimal model. We list the mode identification of the observed periods together with the theoretical periods, 
the stability coefficient cr/, the kinetic energy logE and the rotation coefficient Cu- Also indicated are the relative dispersion of the 
period fit and the amphtude rank of the observed periodicity. 







^ohs 


P,i, 




In? K 




AP/P 


Comments 


/ 


k 


(s) 


(s) 


(rad/s) 


(erg) 




(%) 







3 




112.010 


- 1.070 X 10"' 


40.781 











2 




128.259 


-1.233 X 10-5 


40.751 











1 


155.767 


155.053 


-1.048 X 10"^ 


41.651 




+0.4585 


til 








165.687 


164.995 


-7.363 X 10-^ 


41.741 




+0.4175 


112 


1 


4 




108.498 


-1.404 X 10-5 


40.632 


0.01728 






1 


3 




127.004 


-1.357 X 10^5 


40.710 


0.01367 






1 


2 


149.027 


149.644 


-9.475 X 10"^ 


41.720 


0.03214 


-0.4140 


117 


1 


1 




163.618 


-1.080 X 10"^ 


41.588 


0.01923 






2 


4 




103.956 


-2.518 X 10^5 


40.323 


0.02767 






2 


3 




124.148 


-1.408 X 10^5 


40.690 


0.03422 






2 


2 




138.108 


-2.249 X 10-^ 


41.409 


0.08731 






2 


1 


161.554 


161.986 


-1.390 X 10-^ 


41.487 


0.02621 


-0.2673 


113 


2 





192.551 


192.519 


-8.212 X 10-'° 


44.354 


0.36050 


+0.0168 


114 


3 


4 




100.866 


-3.466 X 10 ' 


40.130 


0.02603 






3 


3 




118.293 


-1.096 X 10-5 


40.785 


0.07479 






3 


2 




130.678 


-7.973 X 10"^ 


40.913 


0.05462 






3 


1 




160.138 


-1.538 X 10"^ 


41.445 


0.04046 






3 







174.449 


-3.289 X 10-* 


42.952 


0.17867 






4 


3 




112.770 


-1.088 X 10^5 


40.767 


0.06969 






4 


2 




128.086 


-1.149 X 10^5 


40.764 


0.03321 






4 


1 


157.581 


157.851 


-1.542 X 10"* 


41.444 


0.05935 


-0.1715 


Its 


4 





168.784 


168.916 


-1.644 X 10-' 


42.302 


0.11349 


-0.0780 


116 


4 


1 




238.583 


-3.504 X 10-'2 


46.073 


-0.01391 







lations also depends on the intrinsic amplitudes, which depend 
on the exact physical conditions inside the star integrated over 
the propagation region of each mode. There is also some ev- 
idence of oscillation power being shifted between modes over 
time, measurable from the relative amplitude changes of the 
observed pulsation periods (insofar as beating between closely 
spaced modes can be excluded). PG 0911+456 is an example 
of this : in the 9 years between the observations of iKoen et al.l 
dHH) and our own, the amplitude of the dominant peak com- 
pared to the second and third periodicity has increased by a fac- 
tor of nearly 3, suggesting a transfer of power into the preferen- 
tial driving of the highest amplitude mode. It remains to be in- 
vestigated whether this is a long-term trend and, indeed, whether 
the amplitude variations reported for other sdB pulsators follow 
any discernable pattern. This will be an interesting incentive for 
future observational projects. 

3.3. Structural parameters of PG 091 1+456 

The optimal model identified leads to a natural determination of 
the three variable input parameters for PG 091 1+456, these be- 
ing the surface gravity log g, the stellar mass M,, and the mass of 
the hydrogen-rich envelope through the quantity log ^(H). In ad- 
dition, we have the effective temperature T^ff from spectroscopy. 
On the basis of these quantities we can derive a set of secondary 
parameters: the stellar radius R (as a function of M* and g), the 
luminosity L (as a function of T^-ff and R), the absolute magni- 
tude My (as a function of g, Teft and M* in conjunction with the 
use of detailed model atmospheres) and the distance from Earth 



d (as a function of apparent magnitude and My). Given that the 
Stromgen magnitude of b -14.6 estimated bv lKoen et al.l(ll999l) 
is associated with large and not readily quantifiable errors, we 
followed the suggestion of the referee and employed the SDSS 
values instead, these being m' = 14.073±0.005, ^'=14.458+0.006 
and r'= 14.885+0.005 (accessible through Vizier@CDS). Using 
the calibration equations given in ISmith et al.l (I2002i) we derive 
a Johnson magnitude of y= 14.663, associated with a formal er- 
ror of 0.007 mags. However, the latter is likely to be severely 
underestimated as the transformation from the SDSS to the 
Johnson-Cousins colour system introduces additional uncertain- 
ties. Comparing the V and B -V values of 4 sdB sta rs cahbrated 
on the basis of the SDSS values to those given by Allard et akl 
(1 19941) reveals typical differences in the magnitudes on the or- 
der of 0.03-0.05 mags. For the relatively faint PG 0911+456 
we therefore adopt y=14.663±0.050, and use this in the compu- 
tation of the distance. Tentative limits could also be set on the 
rotation period Proi and the equatorial rotation velocity Veg (as 
a function of R and Pro/) due to the absence of fine frequency 
structure in the observed period spectrum, assuming that any 
multiplet components would have amplitudes above the detec- 
tion limit. The values and limits for the parameters derived are 
summarised in Table|5] Uncertainties on the three primary aster- 
oseismological quantities log g, M, and log ^(H) were estimated 
from the semi-axes of the 1 cr contours shown in Figures [TOl and 
nn Thes e were calculated acc ordi ng to the recipe describe d in 
detail bv lBrassard et all (1200 Ih and lCharpinet et alj (l2005bl) . 

Comparing the structural parameters inferred for PG 
0911+456 to those determined for other EC 14026 stars from 
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Fig. 12. Comparison of the observed period spectrum of PG 
0911+456 (continuous line segments) with the theoretical pul- 
sation spectrum of the optimal model (dotted line segments) in 
the 100-200 s period range for degree indices ^=0,1, 2,3,4. Note 
that all the theoretical modes plotted are predicted to be excited 
and correspond to low-order p-modes. The radial order k is in- 
dicated below each segment. 

Table 5. Structural parameters of PG 091 1+456. 



Quantity 


Estimated Value 




Teff (K)^ 


31940 


+ 


220 


(0.69%) 




5.111 


+ 


0.002 


(0.03%) 




0.39 


+ 


0.01 


(1.55%) 


log(Me„JM.) 


-4.69 


+ 


0.07 


(1.49%) 


RIRe (M.,g) 


0.134 


+ 


0.002 


(1.00%) 


L/LeiT,s,R) 


16.8 


+ 


1.0 


(4.76%) 


My {g, r,ff, M.) 


4.82 


+ 


0.04 


(0.83%) 


d{V, My) (pc) 


930.3 


+ 


27.4 


(2.9%) 


P,o, (days) 




>68 






Veq (Pro,, R) (knVs) 




<0.1 






' From spectroscopy 



previous asteroseismological analyses, we find that both the total 
mass and the envelope mass are lower than average (see Figure 
[T3]) . In fact, the total mass is the smallest ever measured for 
an sdB star and lies well below the canonical model value of 
~ 0.47 Mq, while s till fall ing int o the low-mass ta il of the dis- 
tribution expected (iDorman etraLlI 19931; iHan et al.l(20()2l |2003|) . 
Assuming that PG 09 1 1 +456 is indeed a single star, the mass de- 
rived lies just below the expected ran ge predicted from the white 
dwarf merger scenario proposed by Han et al.l (12002 ); however, 
the authors specifically state that their values constitute upper 
limits since they assume that all the mass contained in the two 
merged helium white dwarfs goes into the subdwarf. In this con- 
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Fig. 13. A graphical illustration of the total masses and relative 
hydrogen envelope masses determined for the 7 rapidly pulsat- 
ing subdwarf B stars so far submitted to asteroseismology (for 
references see the Introduction). 



text, the authors suggested searching for circumstellar matter 
around sdB's thought to be single stars. A positive detection 
would lower the minimum mass expected for products of the 
WD channel and perhaps explain the small mass determined for 
PG 09 1 1 +45<fl As for the hydrogen shell mass, our results indi- 
cate that it is slightly smaller than that of most EC 14026 stars, 
even when considered relative to the low total stellar mass (see 
Figure [T3]l. This is precisely what would be expected if this star 
did indeed form via the WD merger scenario. 



4. Conclusion 

We obtained 57 hours of broad-band time-series photometry as 
well as high S/N low- and medium resolution time-averaged 
spectroscopy for the EC 14026 pulsator PG 091 1+456. Our ob- 
servations led to refined estimates of the star's atmospheric pa- 
rameters and the detection of 7 independent harmonic oscilla- 
tions, 4 more than were known previously. There was no sign of 
frequency splitting over the 68 -day period during which the pho- 
tometry was obtained, indicating a slow rotation rate. Fixing the 
effective temperature to the spectroscopic value and conducting 
an asteroseismic search in 3-dimensional \ogg - M» - log ^(H) 
parameter space enabled the identification of several families of 



^ Note that the asteroseismic mass determination is robust. We car- 
ried out additional searches in parameter space forcing the mass to the 
canonical value of 0.47 Mq and found that the resulting best fit solu- 
tion presented an unacceptably bad period match with a dispersion of 
< AP/P >=1.5 % or AP=2.5 s (compared to < AP/P >=0.26 % for 
the optimal model). Moreover, the mode identification inferred was less 
than convincing, the only radial mode being assigned to one of the low- 
est amplitude pulsations. The period fit was improved only marginally 
by varying the mass in the 0.46-0.48 Mq range. 
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models that could reproduce the observed periods to within less 
than 1 %. While some of these were rejected from the outset 
due to obvious inconsistencies with the spectroscopic estimate 
of log g or implausible associated mode identifications, we re- 
tained two promising solutions for closer inspection. Unlike in 
some previous asteroseismological studies, it was not immedi- 
ately obvious which model was to be preferred on the basis of the 
structural parameters alone; instead we used the inferred mode 
identification, in particular the degree index {, to discriminate 
between the two. The main difference between the solutions was 
that one identified a relatively high amplitude peak with an i-3 
mode, while the other required only modes with £-0,1,2 and 4. 
Since detailed computations reveal (-3 modes to have extremely 
small disk-integrated amplitudes that would most likely not be 
detectable, we favoured the latter and adopted it as the optimal 
model. 

The inferred structural parameters for PG 09 11 -1-456 include 
the total stellar mass and the thickness of the hydrogen-rich 
shell, two quantities that can normally not be derived using other 
means but are invaluable for a detailed understanding of sub- 
dwarf B stars' evolutionary history. The total mass determined 
is smaller than that found for any EC 14026 star to date and 
places our target at the low-mass end of the predicted distribu- 
tion. Similarly, the hydrogen envelope is measured to be thinner 
than that of most sdB's studied so far If PG 0911+456 is con- 
firmed to be a single star, as is suspected from its slow rotation, 
negligible radial velocity variation, and absence of a compan- 
ion's spectroscopic or near-IR photometric signature, it may be 
the product of a WD merge r accor ding to the evolutionary chan- 
nels proposed by Han et al] (12002 ', 2003). In this case, we would 
indeed expect the hydrogen envelope mass to be smaller than for 
an sdB having undergone a CE or RLOF phase. The low total 
mass derived would tend to support a non-canonical evolution- 
ary history, even if it lies slightly below the mass distribution 
predicted from a WD merger 

An obvious follow-up study for PG 0911+456 is to verify 
the asteroseismological solution found on the basis of additional 
observations. These could aim for a higher S/N level, thus en- 
abling the detection of further pulsations and strengthening the 
constraints on the asteroseismic model. One of the main chal- 
lenges we faced in the search of parameter space was the rela- 
tively large number of models that could account for the oscil- 
lations observed quite accurately; it would be very instructive 
to see whether any newly found periods can also be fit by the 
optimal model isolated. Observations containing wavelength- 
dependent information from which modes may be partially 
identified are another option. Following recent theoretical in- 
vestigations into the amplitude-wavelength dependenc e of a 
mode on its degree ( (Ramachandran et al. (2004), Randall et^ 
(l2005h ). there has been a surge in observational efforts to obtain 
multi-colour photometry of EC 14026 sta rs, most notably us- 
ing the 3-channel CCD ULTRACAM (e.g. lJefferye t al. (2005)) 
and the Mont4kccd predecessor, LaPoune I (e.g. Tbnta ine et al.l 
(12006.) ). While discriminating between low-degree modes with 
{=0,1,2 has proved extremely challenging, {=4 modes exhibit 
a more clearly distinguishable amplitude-wavelength behaviour 
(see e.g. Figure 26 of Randall et al. (2005)). Given the necessary 
data, it should therefore be possible to confirm the identification 
of the two {=4 modes inferred for PG 091 1+456 and thus verify 
the structural parameters computed. 

The work presented here constitutes the 7th detailed aster- 
oseismological analysis of a rapidly pulsating subdwarf B star. 
While we estimate that the structural parameters of around 20 
targets are required to start detailed comparisons with evolution- 



ary theory, first tentative efforts in this direction look promising. 
Nevertheless, it is clear that there is still ample room for im- 
provement on the modelling front. Firstly, the fact that the dis- 
persion between the observed and theoretical periods of the op- 
timal model is generally an order of magnitude higher than the 
measurement accuracy indicates remaining shortcomings in the 
models. We are currently working on full evolutionary (rather 
than envelope) "third-generation" models to address this prob- 
lem. The reliability of the "optimal" models identified during the 
search of parameter space is another issue. Although we do ap- 
ply cross-checks such as compatibility with non-adiabatic theory 
and spectroscopic values of the atmospheric parameters, the lat- 
ter (especially Teff) must often be used to discriminate between, 
or constrain, regions of minimum S ^ and can no longer be em- 
ployed as independent estimates. Moreover, the period ranges of 
unstable oscillations computed from our non-adiabatic pulsation 
code are sensitive mostly to T^^ff and log g and are largely inde- 
pendent of the model mass and envelope thickness. It is therefore 
vital that additional checks are carried out with regard to the ro- 
bustness of the "forward" approach if the structural parameters 
inferred from asteroseismology are to be compared with evolu- 
tionary predictions in a quantitative manner. The most obvious 
way of doing this is by detecting more frequencies from higher 
S/N observations or constraining the identification of the degree 
{ of individual modes from multi-wavelength time-series data. 
Such efforts are ongoing, and will likely prove invaluable for the 
future of sdB star asteroseismology. 
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